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Abstract 

In this work we present a general treatment of a bosonic dissipative network: a chain of coupled 
dissipative harmonic oscillators whichever its topology, i.e., whichever the way the oscillators are 
coupled together, the strenght of their couplings and their natural frequencies. Starting with a 
general more realistic scenario where each oscillator is coupled to its own reservoir, we also discuss 
the case where all the network oscillators are coupled to a common reservoir. We obtain the master 
equation governing the dynamic of the network states and the associated evolution equation of 
the Glauber-Sudarshan P-function. With these instruments we breafly show how to analyse the 
decoherence and the evolution of the linear entropy of general states of the network. We also show 
how to obtain the master equation for the case of distinct reservoirs from that of a common one. 
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I. INTRODUCTION 



Over the last few years interest has grown for better understanding the phenomena of 
coherence and decoherence dynamics in quantum networks, especially in connection to the 
protocols for quantum-state transferring and quantum-state protection for information pro- 
cessing. Beyond the quest for conditions that weaken the system-reservoir coupling , the 
search for mechanisms to bypass decoherence started with quantum-error-correction codes 

and goes through the program of engineering reservoirs |4i]. Moreover, in a closer contact 
with quantum networks, the investigation of collective decoherence resulted on what has 
been called a decoherence- free subspace js], [gI [t]. Interestingly enough, while, in general, 
the protocols for quantum error-correcting codes presuppose that quantum systems deco- 
here independently, the decoherence-free subspace is generated by distinct quantum systems 
coupled to a common reservoir. 

Regarding state transfer, the controlled coherent transport with splitting of atomic wave 
packets Sj, and the evolution of macroscopically entangled state [Sj have been analyzed 
within the context of optical lattices. The dynamics of Bose-Einstein condensates in a one- 
dimensional optical lattice is also investigated and a class of spin networks have been 
proposed for perfect state transfer of any quantum state in a fixed period of time ll|. In 
the context of systems of coupled harmonic oscillators (HOs), which we focus in the present 
work, the dynamics and manipulation of entanglement was analyzed in Ref. [12|]. Evidently, 
the pressure for the implementation of logical operations with an increasingly larger number 
of quantum systems will decisively drive on the quest for controlled coherent transport in 
quantum networks. 

Concerned with a simple network of two coupled resonators, Raimond et al. 13| have 



presented the blueprint of an experiment in which the decoherence of a mesoscopic superpo- 
sition of radiation states becomes a reversible process. A theoretical model of the proposal 



in Ref. |13|] is given in Ref. 14|, where the coupling of the resonators to their environ- 



ments is taken into account when the reversibility of coherence loss is analyzed. In Ref. 
[isl . the authors assume that only one of the resonators in Refs. [isl, Q] is interacting 
with a reservoir to derive a master equation in the case where the resonators are strongly 
coupled. It is shown that the relaxation term is not simply the standard one, obtained by 
neglecting the interaction between the cavities, i.e., dissipation is not additive for strongly 
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coupled systems. Finally, in Ref. [l6| both resonators are considered to be lossy, as in Refs. 



13 



14| . and the regime of strongly coupled cavities is also analyzed, as done in Ref. [15(]. A 



detailed analyses of the coherence and decoherence dynamics of quantum states is presented 



m 



161], regarding this network of two coupled resonators, including a study of the correla- 



tion between the fields in both resonators through the excess entropy. The phenomena of 
electromagnetically induced transparency and dynamical Stark effect are also analyzed in a 
network of two coupled dissipative resonators jiyj ]. 

In the context of complex networks, composed by a large number of subsystems, in 



Refs. Q 



the authors present a detailed treatment of the coherence and decoherence 



181 ] a symmetric network of 



dynamics in arrays of coupled dissipative resonators. In Ref. 
interacting resonators are considered, where each oscillator interacts with each other, apart 
from its own reservoir. A different topology is analyzed in Ref. [3], where a central oscillator 
is assumed to interact with the remaining A^ — 1 peripheral and noninteracting oscillators. In 
both topologies, the decoherence process is analyzed by focusing on a single resonator which, 
apart from interacting with its own reservoir, also interacts with the remaining A^ — 1 coupled 
resonators plus their respective reservoirs. Considering all resonators with the same natural 
frequency t^o and all couplings with the same strength A, master equations are derived 
for both weak (A ^ ujq) and strong (A ~ uq) coupling regimes. From such development, a 
detailed analyzes of the emergence of relaxation- and decoherence-free subspaces in networks 
of weakly and strongly coupled resonators is presented in Ref. . The main result in Ref. 



201 ] is that both subspaces are generated when all the resonators couple with the same 



group of reservoir modes, thus building up a correlation (among these modes), which has 
the potential to shield particular network states against relaxation and/or decoherence. 

It is worth noting that recent results regarding entanglement and nonclassical effect in 
collective two-atom systems 2l|], retains some resemblance with those discussed above for 
networks of coupled resonators. Beyond the entanglement dynamics which is a crucial but 
recurrent ingredient of any network, the collective damping effects coming from two-atom 
systeni s I21I] c an be directly identified with those in a network of dissipative oscillators 



3,Q, 



19|, 



201 ] . Such collective damping effect are certainly in the basis of the nonadditivity 



of decoherence rates observed in the network of dissipative oscillators Q, 1^, 1^, 
as in superconduc ting qu bits |22|. 



20| as well 



Since in Refs. 



18 



191 ] two different topologies are analyzed independently, for the par- 
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ticular case where all resonators have the same natural frequency and all couplings have 
the same strength A, in this work we present a unified approach for treating a bosonic dis- 
sipative network. Such approach holds for whichever the topology of the network, i.e., for 
whichever i) the way the resonators are coupled among them, ii) their coupling strengths 
and in) natural frequencies. 

In the Section II, towards the derivation of the master equation governing the evolu- 
tion of the network, we present our model. Considering first a nondissipative network, we 
show how to derive particular topologies from the general case of a symmetric network 
where each oscillator interacts with each other. In Section III, the evolution equation of 
the Glauber-Sudarshan P-function is obtained as a c-number map of the master equation 
in operator form. Thus, in the context of dissipativc networks, we show how to derive 
particular topologies from a general symmetric dissipative network where each oscillator is 
coupled to its respective reservoir apart from interacting with each other. Solutions in terms 
of the Glauber-Sudarshan P-function, for general initial states of the network, are given in 
Section IV together with a brief analysis of decoherence and the hnear entropy. Finally, the 
concluding remarks are presented in Section VI. 

II. GENERAL TREATMENT OF A BOSONIC NETWORK 
A. The model 

Setting from here on that the indexes m, m', n, and n'run from 1 to A^, we start considering 
a general Hamiltonian for a bosonic network, H — Hs -\- Hr -\- Hr, involving a network of N 
coupled oscillators 




(1) 



N distinct reservoirs, composed by a set of /c = 1, . . . , oo modes. 




(2) 



m 



k 



and the coupling between the network oscillators and their respective reservoirs 




(3) 



m 



k 
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^Ink i^mk) IS the creatioii (annihilation) operator for the kth bath mode cumk coupled to 
the mth network oscillator Um whose creation (annihilation) operator reads aj^^ {am ). The 
coupling strengths between the oscillators are given by the set {Xmn}, while those between 
the oscillators and their reservoirs by {Knfc}- Before addressing the dissipative process 
through Hamiltonians ([2]) and ([3]), we focus first on Hamiltonian Hs to show how to derive 
different topologies of a nondissipative network of coupled harmonic oscillators. Rewriting 
Hs in a matrix form Hs = hJ2m n 'An^mnO'ni its elements are given by 

{Um '^^ rn = n 
1 (4) 

whose values characterize whichever the network topology: the way the oscillators are cou- 
pled together, the set of coupling strengths {Xmn}, and their natural frequencies {wm}- 



B. From the general matrix H to particular nondissipative topologies 



To illustrate the procedure to construct particular nondissipative topologies we consider 
four different cases: the i) symmetric, ii) central, in) circular, and iv) linear networks. For 
the case of a i) symmetric [sym) network, sketched in Fig. 1(a), all the oscillators are 
coupled together, with all matrix elements of Ti being not null 

f LOi Ai2 Ai3 ■ ■ ■ AiAT \ 
Ai2 UJ2 A23 ■ ■ ■ A2Ar 
Al3 A23 ^3 ■ ■ ■ AsAT 



n 



sym 



(5) 



\ Aiiv A2Ar X3N ■ ■ ■ I 

In a ii) central (cent) network, sketched in Fig. 1(b), only one selected oscillator, the central 
one, is assumed to interact with the remaining — 1 noninteracting peripheral oscillators. 
Labeling the central oscillator by 1, with the peripherals running from 2 to N, the matrix 
Ti. has the form 

f UJi A12 A13 ■ ■ ■ AiAT \ 

A12 CJ2 ■ ■ ■ 

Ai3 CU3 ■■■ , (6) 



n 



cent 



IN 







UN j 
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where only in first column and row all the nondiagonal elements are not null. As depicted 
in Fig. 1 (c), in a iii) circular {circ) network the kth oscillator is coupled to the {k ± l)th 
oscillators, with the additional condition that the iVth oscillator be coupled to the first one. 
The matrix 7i is given by 

/ Ui Ai2 ■ ■ ■ AiAT \ 
Ai2 UJ2 A23 ■ ■ ■ 

X23 UJ3 ■■■ 



He 



(7) 



\ AiAT ■■■ UN / 

Finally, the iv) linear (lin) network follows directly from the circular one by turning off the 
coupling between the first and the A^th oscillators. The matrix Ti. obtained for this case has 
the three-diagonal form 

I UJl \l2 ••■ \ 
A12 1-1^2 A23 • ■ ■ 

A23 CU3 ■ ■ ■ 



n 



lin 



(8) 



\ ■■■ tujv/ 
Next, we treat the general situation of the dissipative network, described by Hamiltonian 
H = Hs + Hfi+Hj, through the standard perturbative approach in the system-bath coupling 
strengths {Knfc}- Since all the particular topologies follow from the case of a symmetric 
network, choosing appropriately the elements of Ti, we shall obtain the reduced density 
operator of the coupled oscillators from this general topology. 



C. The master equation of a bosonic dissipative network — Direct and indirect 
dissipative channels 



To obtain the master equation of the network we first diagonalize the Hamiltonian 7i 
through a canonical transformation 



(9) 



where the coefficients of the mth line of matrix C define the eigenvectors associated to 
the eigenvalues of matrix ([5]). With C being an orthogonal matrix, in that = 



C ^, it follows the commutation relations [Am,^.]^] = Smn and [y4m,^n] = 0, enabling the 
Hamiltonian H to be rewritten as H = Hq + V , with = AnCnm and 



m m k 

y = ^J2Y1 '^nmVmk (^Lfc^n + &mfc^l) • (10b) 
m,n k 

With the diagonalized Hamiltonian Hq we are ready to introduce the interaction picture, 
defined by the transformation U{t) = exp {—iHot/h), in which 

V{t) = hJ2 {Omnmi + OUt)An) , (H) 

m,n 

where Omn(t) = CnmJ2k'^'m.k^w[~'i i^mk ~ ^n)t]^rnk- Ncxt, wc assumc the interactions 
between the resonators and the reservoirs to be weak enough in order to perform a second- 
order perturbation approximation followed by tracing out the reservoir degrees of freedom. 
We also assume a Markovian reservoir such that the density operator of the global system 
can be factorized as pi...Ar(t)®p/?(0). Under these assumptions we obtain the reduced density 
operator of the network of dissipative coupled resonators given by 



fdt'TTR[V{t), [V{t'),pn{0)®p,_N{t)]] . (12) 
Jo 



dpi,...,jv(t) ^ 1_ 

dt ~ h? 

Since for a thermal reservoir {bmkbnk') = (^InkKik'^ ~ 0' have to solve the integrals 
appearing in Eq. (fT2l) . related to correlation functions of the form 

/ dt' (^Ol^n{t)Om'n'{t')) = CnmCn'm' / d t' VmkVm'k' {b^mk^rn'k' / 

X exp {i [{Umk - ^n) t - {uJm'k' " ^n') t']} ■ (13) 

Considering that the reservoir frequencies are very closely spaced to allow a continuum 
summation and defining the average excitation of the mth mode associated to the kth 
reservoir (i^) as (bl^ki^)bm'k'{j^')^ = 27r5„m'N„(z/)5 {u - u'), we obtain 

' dt' {Ol^{t)Om'n'{t')) = Smm'CnmCn'me'^^-'-^-^' 

/■* f°° d V 

X dt'y^ _[Vm{y)am{y)f^m{y)e-'^''-^-'^^''-'\ (14) 
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with am{i^mk) being the density of states of the mth reservoir. Assuming, as usual, that 
V^(ron), (Tmi'cOr,) and N^(ron) are slowly varying functions, we obtain after the variable 
transformations e — v — w^i and t — t — t' the simplified form 



'■* N 



J ^t' {^Ln{'t)^m'n'{t')) = '-^^mm'CnmCn'm'ym{'^n')'^m{'^n') exp [i {Wn' - Wn) t] , (15) 

where we have defined the damping rates as 

1 2 /"°° 

Imi'^n) ^ j;^\ym{-^n)(7m{-^n)] I 5{e)de. (16) 

J —Rn 

Back to the Schrodinger picture and to the original field operators a^, we finally obtain 
from the steps outlined above, the master equation 

— [Pl,...,N[t), -tio] + Y <-^n'm(-^nn'7m(^n'j 

m,n,n' 

X {N^(tx7„,) ([a|,pi,...,iv(t),a^] + [aj„, Pi,...,iv(t)an]) 

+ {Nm{rUn') + 1) {[anPi,...,N{t),al^] + [um, Pi,...,Nit)al])} . (17) 



Prom here on we shall focus on the case of reservoirs at OK, leaving for the last but one 
section a brief analysis of the effect of finite temperatures. Defining the effective damping 
matrix whose elements are 

{'UJn')Cn'n, (18) 

n' 

the master equation for the reservoirs at OK simplifies to the generalized Lindblad form 

^J2^{[anPi,...,Nit),al] + [am, Pi,...,Nit)al]) 

m,n 

' {pi,...,N{t)^Ho]+J2^mnPl,...,N{t)^ (19) 



h 



m.n 



where C„mPi iv(^) are the Liouville operators accounting for the direct (m = n) and indirect 

(m 7^ n) dissipative channels, respectively. Through the direct dissipative channels the 
oscillators lose excitation to their own reservoirs, whereas through the indirect channels 
they lose excitation to all the other reservoirs but not to their own. We observe that for 
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Markovian white noise reservoirs, where the spectral densities of the reservoirs are invariant 
over translation in frequency space, such that 7m(^^n') = 7m, the expression f|T8l) reduces to 
Tmn = N'^m^rnn- For thcsc particular reservoirs the indirect channels disappears. We also 



observe that, in the weak coupling regime where {Xmn} ^ {^m'} [18|, ll9j and consequently 
7m(ro„/) ^ •y^^Ujn), we obtain F^n = N'jrn{^m)Smn, such that the indirect channels again 
disappears. Therefore, it is worth nothing that the indirect channels play a significant role 
only in the strong coupling regime where {Xmn} ~ Wm'}- 



III. THE GLAUBER-SUDARSHAN P FUNCTION 

The evolution equation for the Glauber-Sudarshan P function, derived from the master 
equation (fT9l) . is given by 

= 2^ I ^ + l^^rnnVn-Q— + C.C. 1 Pl,...,N{{Vm'} , t) , (20) 

m \ n / 

where we have defined the matrix 7^^, with the elements 

^mn = ^mn/'^ + ^^mn, (21) 

thus generalizing the former matrix Ti ^ to account for the dissipative {D) process. With 
the transformation -Pi,...,Ar({?7m'}, i) = -Pi,...,Af({^m'}, i) exp (^^Fmrn^), and assuming a solu- 
tion of Eq. (!20|) of the form P({?7„} ,t) = P{{rin(t)}), we obtain the differential equation 

dPl,...,N{{Vm'{t)}) _ V- / d7]m{t) 9 \~ ^ r„ 

= J2 ij2'^^nVnJ^ + C.C.] Pi_N{{Vm'},t), (22) 
m \ n / 

which makes possible to calculate the time evolution of the parameters rjmit) through the 
physical quantities of the system appearing on the elements as 

^ = (23) 

n 

Through the transformation rjm(t) = J2n^mnVmit), we diagonalize the matrix thus 
reducing Eq. ( l23l) to the diagonal form drfm{t)/dt = QrnVm, whose solution is rjmit) = 
^mGxp {Qmt)- Therefore, back to the parameters rim{t) we obtain 



m 



(t) = DmnVnit) = ^ Dmn exp (^]„^) An, (24) 



where the elements of the mth column of matrix D define the mth eigenvector associated 
to the eigenvalue Qm of matrix H^, and by setting the initial condition rjmit = 0) = r/^, we 
verify from Eq. that J2n Dmn^n = Vm- Therefore, An = D~^t]^ and, consequently 

r/„(t) = J2 Dmn exp D-^i^l,, (25) 

m' ,n 

leading to the solution for the Glauber-Sudarshan P-function 

Pl,...A{Vn]^) = exp (^T^rnt^ Pl,...,N{{Vn},t = 0) 1 (j)} • (26) 

Therefore, having the P-function at time t = 0, we immediately obtain it at any other time 
by substituting the set {?7„} by {Tin{t)}. 



A. Prom the general matrix Ti^ to particular dissipative topologies 



In subsection II. B we illustrate how to construct particular nondissipative topologies 
from a general symmetric network described by matrix ([5]). Now, after introducing the 
generalized matrix Ti.^ we are in the position to enlarge the focus by constructing networks 
entirely composed of dissipative oscillators or, in a more general fashion, composed by mixed 
nondissipative and dissipative oscillators. Back to the symmetric, central, circular and linear 
networks, when considering that they are all composed by dissipative oscillators, each one 
coupled to its respective reservoir, we obtain for the matrices Ti.^ exactly the same structure 
as those in subsection 2.B. However, from Eq. fl2T]) we verify that the matrix elements Ti^n 
follow from those of Timn multiplied by the imaginary i apart from the correction Fmn/S 
coming from the dissipative process. As an example, for the symmetric network composed 
entirely by dissipative oscillators, as sketched in Fig. 2 (a), the Hamiltonian assumes 
the form 

/ Til Fis ■ ■ ■ FiTv \ 

F21 F22 F23 ■ ■ ■ F2Ar 
F31 F32 F33 ■ ■ ■ F37V 



njD 



^T~^sym ~l~ 2 



(27) 



\ Fati F7V2 FAr3 ■ ■ ■ Tnn ) 

Let us consider, instead, the case of a mixed symmetric {mix — sym) network, composed by 
an even total number of oscillators N, where those designated by odd (even) numbers are 
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nondissipative (dissipative), as sketched in Fig. 2 (b). In this case, the Hamiltonian Ti.^ is 
given by 

/ ••• \ 

T22 ■ ■ ■ r2Ar 
••• 



n 



D 

mix—sym 



in 



1 



sym + 2 



(28) 



\ Tati ^N3 ■ ■ ■ '^NN ) 

As a last example of a mixed network we consider the topology sketched in Fig. 2 (c), 
where three linear chains of coupled oscillators are connected together, with the central 
(lateral) chain(s) being composed by nondissipative (dissipative) oscillators. In this case, 
labeling the oscillators as in Fig. 2 (c), we obtain 

\ 



n 



D 



\ 



( 



+ 



fii Tis 



(29) 



I UJx A12 Ai4 

A12 l-o'2 A23 A25 

A23 cus A36 

Ai4 a;4 A45 A47 

A25 A45 (^5 A56 A58 

-^36 A56 cue Agg 

A47 co'7 A78 

A58 A78 1x^8 -^89 

A89 CUg j 

Therefore, for a given topology, the diagonalization of the matrix HP ^ leading to the 
eigenvalues fi^ and eigenvectors composing the matrix D, enables us to obtain the evolved 
Glauber-Sudarshan P-function and, consequently, the reduced density operator of the net- 
work, given by 



r2i r22 r23 



r4i r42 r43 
000 

Tei Tes 



Tsi r82 r83 

y Tgi r92 Fgs 



19 

r29 



r4g 


r69 



r89 

Fqq 



Pi,...,iv(^) 



(30) 



IV. INITIAL STATES, DECOHERENCE, AND LINEAR ENTROPY 



Starting from two different initial states of the network, given by general superpositions 
of coherent and Fock states, in this section we obtain the evolved reduced density operator 
fl30|) for both cases. 
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A. A general superposition of coherent states 



Let us assume that the initial pure state of the network is given by 

Q 

p{o)=^^'J2^rKm})m}\, (31) 

r,s=l 

where N" is the normahzation factor, is the probabihty amphtudes of the product state 
|{/^m}) — 'S)m=i \Prn)y labcls r and s run from 1 to the integer Q. The superscript r 

stands for the rth state of the superposition while the subscript m stands for the coherent 
state of mth oscillator. We stress that the discrete sum of product states in Eq. f l3ip 
can be substituted by the continuum sum |V'(0)) = A/" / d9A {6) \ {(3m {d)}) with no further 
complication. After some algebra, we verify that the state fj3T|) evolves to the A^-oscillator 
density operator 

where the excitation of the mth oscillator is given by 

Cmit) =Y.^Ut)f^l (33) 

n 

with the time-dependent matrix elements 

Qmn{t) = ^ Dram' ©Xp (34) 
m' 

For the reduced density operator of the mth oscillator we obtain 

Pm{t) = E ^^'^^^^y^y , (35) 

where the influence of all the other oscillators of the network is present explicitly in the 
product {{Pn} \{Pn}) implicitly in the states |Cm(^))- 



B. A general superposition of Fock states 

In the Fock basis we consider the state 

l<^(0)) = 5Z '^ni,...,n„\nu...,nN) , (36) 

ni,...,njv 
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where rim stands for the photon number of mth oscillator and Cni,...,nN is the probability 
amplitude of each state in the superposition. After a lengthy calculation we verify that the 
iV-oscillator density operator is given by 

IIJL (r, -a Vk r 

m qm=0 fcm=0 

X \J^i{qi} , {ke} , t)) {T{{me -n, + q,} , {h} , t) \ (37) 

where 



Um + km)\ '"'^'^ {QmAt)y-^'--^-'' (e„,;v(t)) 



f-m.N-1 



'''' Jm — ^ /^m,l — u 



r=l 



where 5 (x) equals unity for x = 0, being null otherwise. 



C. State transfer and recurrence dynamics 

From the reduced density operators pmit) following from Eqs. ( l32il and ( l37|l it is directly 
to verify the transfer of an initial state prepared in the mth oscillator to the remaining one 
of the network, followed by the recurrence of this state back to the mth oscillator. The 
probability of recurrence of an initial state Pm(0) prepared in the mth oscillator is given by 
the expression 

Vnit) = Tr [pm{t)Pmm , (39) 

which is also a measurement of the fidelity of the initial state Pm,(0), expected to decrease 
due to the dissipative process. For the probability of transfer of the initial state Pm(0) to a 
particular nth oscillator picked up from the remaining — 1 of the network, we get 

VT{t) = Tr [pnit)pmm . (40) 



From Eqs. ( I39i) and ( HOl) it can be verified — as analyzed in details in Refs. ISj and 
[lol for the particular symmetric and central topologies, respectively — that an initial su- 
perposition prepared in the mth oscillator bounces between its original oscillators and all 

13 



the remaining oscillators of the network. Evidently, the dynamics of a given prepared state 
through the network can be manipulated through the choice of the topology. 



D. Decoherence 

From a given initial superposition of coherent states (l3T|l and the density operator of 
the network fl5^ . we can estimate the decoherence time of an arbitrary chosen off-diagonal 
element of the density operator relatively to the relaxation time of the diagonal elements. In 
fact, the literature concerned with the decoherence of A^-dimensional superpositions acquaint 
only for relative time-decay measurements by which the larger the distance from the main 



diagonal of the matrix elements, the smaller are their decay time [2J]. From Eq. (!32|) 
we verify that such relative time-decay measurements follows from the real part of the 
coefficients / ({Cm('^)} l{Cm(^)}) which is directly computed from the initial 

state of the network together with Eq. ( I33l) . 

However, additional ingredients concerning the decoherence dynamics arise when consid- 
ering a network of dissipative quantum systems. In Ref. [l^, where a minimal network 
of two dissipative oscillators is analyzed, it is demonstrated that the decoherence time of 
a "Schrodinger cat"-like state prepared in one of the oscillators can be doubled compared 
to that when the same state is prepared in an isolated dissipative oscillator. This result 
follows when the decay rate of the oscillator, where the state is prepared, is significantly 
larger than the other one composing the network. A generalized analysis of decoherence for 
the case of a symmetric network of dissipative oscillators is presented in Ref. 20j], where 
the physical ingredients that enable the emergence of relaxation-free and decoherence-free 
subspaces are exposed. On this regard, a detailed study of the optimum topologies leading 
to maximum decoherence times of superposition states prepared in particular oscillators of 



dissipative networks will be presented elsewhere 231]. The memory devices presented in Ref. 



231], which follows from the general formalism presented here, combines both ingredients: i) 
the large decay rate of the storage oscillators of the network — those except the one where 
the state to be protected is prepared — and ii) specific dynamics of this state through the 
network, achieved by engineering particular topologies. 
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E. Linear entropy 



From the density operator in Eq. (132|) . we are able to calculate the linear entropies for 
the mixed states of the whole network, 5i^...^7v(^); of oscillator 1 (or any other particular 
oscillator), Si{t), and of all the remaining — 1 oscillators, S2,...,N(t). which are given by 

<Sl,...,Nit) = 1 - Tri,...,Ar [pl,...,Nit)f 



r,s,p,q 



X exp 



(41a) 



5i(t) = l-Tri [p^{t)f 



r,s,p,q 



X exp 



(41b) 



<S2,...,N(,t) = 1 - Tr2,...,Ar [p2,...,N{t)f 

N 



r,s,p,q 



X exp 



(41c) 



where the second equality follows from the initial state flHTl) and, as like as r and s, the labels 



p and q run from 1 to the integer Q. With these expressions we can analyze, as discussed 
m Refs. [l8|,Q,Q, 

the evolution of the correlation between the reduced state of oscillator 
1 and that of all the remaining — 1 oscillators, through the excess entropy, defined as 



Sit) ^ Si{t) + S2,...,N{t) - 



(42) 



The excess entropy can also reveals (through a residual value of E{t)) the development of 
an inevitable background correlation between all the network oscillators which thus become 
permanently entangled 

[l8|,Q,Q- 

This background correlation arises from two different 
mechanisms: First, the action of the indirect channels {LmnPi,...,N{t)) and, secondly, the 
action of the direct channels {LmmPi,...,N(t)), when the decay rates Tmn are different from 
each other. For equal decay rates, the direct decay channels do not contribute to the 
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development of the background correlation. The indirect channels thus play an important 
role in the entanglement process in dissipative networks. 



V. THE CASE OF A COMMON RESERVOIR FOR THE WHOLE NETWORK 

In this section we extend our analysis to contemplate the case where all the oscillators 
of the network are coupled to one common reservoir (at OK). For a particular symmetric 
network, the different results following from both cases of distinct or a common reservoir 
have been discussed in Ref. 20|] in connection with the emergence of relaxation-free and/or 



decoherence-free subspaces. In Refs. [16|, uM, llSl, |20| a brief discussion is also provided about 
the rather unusual scenario of a common reservoir. 

The Hamiltonian for the case of a common reservoir is given by 

+ ^ uJkblbk + ^ ^ ^ Knfc (blam + haln^ ■ (43) 

k m k 

Following the same steps as in Section II, we obtain from Hamiltonian fH5]) the same master 
equation (IT^ derived previously for the case of distinct reservoirs, but with the effective 
damping matrix (ITSl) replaced by 

^mn = N imm'{'^n')Cn'm'Cn'n-, (44) 

m',n' 

where 

ft 1'°° (\v 

Un{wm')= / dr / —a\v)Vm{v)Vn{v)e-'^^-^^'^^'-^\ (45) 



^0 



TT 



The correlation factor S,mn arises from the fact that both network oscillators, m and n, 
may interact indirectly through their common reservoir. To analyze more closely such a 
correlation, we assume (as usual for the case of weak coupling between the system and 
the reservoir), that the network oscillators only interact with the reservoir modes in the 
neighborhood of their normal modes. Under this assumption the maximum correlation 
takes place when both oscillators m and n are identically coupled with the same group of 
reservoir modes, i.e., when Vm{v) = Ki(^)- Otherwise, a partial correlation arises when the 
coupling between the oscillators with the reservoir modes turns out not to be identical, i.e., 
{Kn('^)} n {Vn{i')} 7^ 0. In this case, the oscillators may still be coupled with the same 
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group of reservoir modes, but with different strengths, or be coupled with different groups 
of reservoir modes apart from a common intersection of them. The correlation between 
the oscillators disappears only when {Vm{i^)} H = 0, i-e., there is practically no 

intersection of common reservoir modes coupled to both oscillators. 

It is particularly interesting to note that in the case where {KTi(^)} H = and 

consequently ^mn = 0, only the self-correlation survives, which reduces to the damping 
factor 6mn1m in Eq- f|T6|) . apart from the unique frequency distribution o"(z/) of the common 
reservoir. By its turn, when assuming the coupling strengths Vm between the oscillators 
and the common reservoir to be all different to compensate the unique o"(z/), the effective 
damping factor arising from the self-correlations reduces to that of the case of distinct 
reservoirs in Eq. flTS]) . Therefore, it is possible to derive the master equation for the case of 
distinct reservoir starting from that of a common one, under the condition that no correlation 
between two oscillators is induced by their common reservoir. Conversely, it is also possible 
to shift from the case of distinct reservoir to that of a common one assuming that all the 
reservoirs presents the same frequency distribution cr(z/) and the limit of strong interactions 

nrrnn 

between the network oscillators. In this limit, as discussed in Refs. |16l . 1181 . I19I . |20|, the 
condition NX^n ^ ^^m' niust be satisfied. The interesting aspect of such a condition is that 
it can be fulfilled for coupling strengths Xmn ^ t^m', as long as a sufficiently large network 
is provided (A^ > uJm'/Xmn)- 



VI. THE CASE OF RESERVOIRS AT FINITE TEMPERATURES 

Since a formal approach for the case where the reservoirs are at finite temperatures is 
somewhat demanding, here we shall present a brief qualitative analysis of this case. Our 
analysis focus on the normal-mode oscillators Wm (represented by the operators and 
Al^) under the assumption that all the coupling strengths between the original oscillators 
ujm (represented by am and a}^) and their respective reservoirs are around the same. We 
first note that the normal-mode oscillators are decoupled from each other, whereas each 
one interacts with all the N reservoirs as described by Hamiltonian in Eq. ffTOj) . Evidently, 
when the couplings Xmn between the original oscillators are all turned off, the normal-mode 
oscillators degenerate into the original one. Moreover, when the coupling strengths are 
significantly smaller than the natural frequencies of the original oscillators, i.e., N {Xmn} ^ 
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{cUm'}, the magnitude of the interaction between the mth normal-mode oscillator with the 
mth reservoir is significantly larger than those with the remaining N — 1 reservoirs. In 
fact, when N {X^n} {^m'}, the indirect channels are not quite effective as pointed out 
above. Otherwise, when N {\„in} ~ {^^m'}? the indirect channels become as effective as the 
direct one, and the magnitude of the interactions between the mth normal-mode oscillator 
with all the reservoirs becomes quite the same. From the above qualitative observations we 
next discuss the effects of temperature in both cases of distinct reservoirs and a common 
reservoir. 

A. Distinct reservoirs 

For the case of N distinct reservoirs at finite temperatures Tm, in the regime where 
N {Xmn} -C {uim'}, wc thus concludc that the (non mill) mean energy (E^) of the mth 
normal-mode oscillator will practically be defined by its associated mth reservoir. Conse- 
quently, in the steady state configuration, each of the normal-mode oscillators presents a 
different mean energy which is defined by the equilibrium reached with all the reservoirs, but 
mostly with its associated reservoir. As expected, in the regime where N {Xmn} ~ {^m'}, all 
the normal-mode oscillators present approximately the same mean energy, since the magni- 
tudes of their couplings with the difi'crcnt reservoirs are approximately the same. Summariz- 
ing, in spite of the different temperatures of the reservoirs, in the steady state configuration 
of the regime where N {Xmn} ~ {<^m'}, all the normal- mode oscillators vibrate with approx- 
imately the same mean energy, whereas in the regime N {Xmn} <S {^m'}, each normal-mode 
oscillator exhibits a different mean energy. 

B. A common reservoir 

When a common reservoir is considered, the scenario is much like that arising in the case 
of distinct reservoirs in the regime N {Xmn} ~ {^m'}- In fact, as the coupling strengths 
between the normal-mode oscillators and their common reservoir are around the same, in 
the steady state configuration all the normal-mode oscillators exhibt around the same mean 
energy. 
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VII. CONCLUDING REMARKS 



Motivated by the necessity to better understanding the coherence and decoherence dy- 
namics of quantum states in networks composed by a large number of dissipative quantum 
systems, we have studied chains of dissipative harmonic oscillators. We presented previous 
results related to different topologies, by starting with the simplest case of only two cou- 
pled oscillators [IGj] and generalizing the analysis to the case of coupled oscillators in a 



symmetric 



18] and a central-oscillator network 19]. The developments in Refs. 



were extended to the anal ysis of the physical ingredients responsible for the emergence of 
decoherence-free subspace [20]. 



Since in Refs. 



18 



19| we have treated two particular topologies independently, in the 



present contribution we have presented a general formalism to treat whichever the topology 
of a chain of dissipative harmonic oscillators. Starting from a symmetric network, were all 
the oscillators are coupled together, apart from being coupled to their respective reservoirs 
(or to a common one), we have derived the master equation and the associated evolution 
equation of the Glauber-Sudarshan P-function. We thus have showed how to particularize 
such results for whichever the specific network. We also shown how to obtain the master 
equation for the case where each oscillator is coupled to its respective reservoir starting from 
that where all the oscillators are coupled to a common reservoir. 

The presented formalism is quite general and can be used to compute the decoherence 
time of pure or mixed states prepared in a particular oscillator of the network or even in a 
cluster of oscillators of the network. The correlation between the states of the network can 
also be computed through the excess entropy defined for a bipartite system 16, [l8|, [1 



It is worth stressing that an exact model to treat a dissipative quantum oscillator can be 
extracted from the present treatment. To this end, we have to pick up a single oscillator 
of the network — our system of interest — and couple it to all the other oscillators which 
play the role of the reservoir. Thus, disconsidering the reservoirs which we have treated here 
through the standard perturbative approach, we end up with an exact formalism to account 
for the dissipative effects over a harmonic oscillator. Evidently, such an exact formalism 
can also be pursued for dissipative systems other than a harmonic oscillator. An extensive 
analysis of the effects coming from the exact treatment of dissipation emerging from this 
work will be presented elsewhere. 
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Figure captions 

Fig. 1 Sketches of nondissipative symmetric (a), central (b), circular (c), and linear (d) 
networks. 

Fig. 2 Sketches of a dissipativc symmetric network (a), a mixed-symmetric network (b), 
composed by dissipativc and nondissipative oscillators, and a mixed network composed of 
dissipative and nondissipative chains of oscillators (c). 



21 



Fig. 1(a) 




Fig. 1(b) 




Fig. 1(c) 




Fig. 1(d) 




Fig. 2(a) 




Fig. 2(b) 



Fig. 2(c) 




